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1. Introduction 

The Kaon B parameter 



(Tf\ s -Y^l-Y5)d-sU\-Y5)d\K ) 
(8/3) (1° |^75 d\0) (0\sYpKd\K°) 



is one of the fundamental weak matrix elements which have to be determined theoretically for de- 
ducing CP violation phase of the Cabibbo-Kobayashi-Maskawa matrix from experiments. Lattice 
QCD calculation may be an ideal tool to determine the matrix element precisely from the first prin- 
ciple. An essential step towards precise determination of B K is to control systematic error in the 
renormalization factor. Recently non-perturbative renormalization factor is preferably employed to 
remove errors in the perturbative one. Among several non-perturbative schemes on the lattice the 
Schrodinger functional (SF) scheme [|l|] has an advantage that systematic errors can be evaluated 
in a controlled manner. 

A few years ago the CP-PACS collaboration calculated Bk using quenched domain-wall QCD 
(DWQCD) with Iwasaki gauge action [Jsj] . Their result was renormalized perturbatively at one loop 
and have shown a good scaling behavior with small statistical errors. A main purpose of this paper 
is to derive a non-perturbative renormalization factor 3fg K which convert the bare Bk of the CP- 
PACS collaboration to the renormalization group invariant (RGI) Bk- We adopt the SF scheme 
as an intermediate scheme to avoid systematic uncertainties due to the finite lattice spacing. The 
renormalization factor 3fs K (go) is given at a fixed bare coupling and its non-perturbative evaluation 
is decomposed into three steps in the SF scheme as 

%B K {go) =ZyJ +AV (oo,/X max )Zyj > +Ay (/X max ,/X m i n )ZgJ(go,Mmin)- (1-2) 

We start from a renormalization factor at a low energy hadronic scale /i m i n 

where Zw+aa is a renormalization factor for the parity even part of the left-left four-quark operator 
and Za is that for the axial vector current. A reason why we define the renormalization factor at 
the hadronic scale is to suppress lattice artifacts by a condition a/Vm *C 1. This factor depends on 
both renormalization scheme and lattice regularization. Multiplying it with a lattice bare operator, 
the regularization dependence is canceled and only the scheme dependence remains. 

Z\/a+av (/Vax>/Vix) represents non-perturbative RG running of the parity odd part of the left- 
left four-quark operators from the low energy scale /i m j n to a high energy scale /i max = 2 7 jU min , 
where perturbation theory can be safely applied. Among three steps this requires the most exten- 
sive calculation. Since this factor evaluated in the continuum limit does not depend on a specific 
lattice regularization, we can employ ZyJ +Ay (/i max ,/i m i n ) evaluated previously by the Alpha col- 
laboration with the improved Wilson fermion action[||], instead of calculating it by ourselves. Even 
though renormalization factors for the parity even and parity odd parts differ on the lattice without 
chiral symmetry, the difference disappears in the continuum limit, where the chiral symmetry is 
recovered. 
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The last factor Z^ A+AV (oo, jU max ) is the RG evolution from the high energy scale /x max to infinity, 
which absorb the scale dependence to give the RGI operator. Since we are already deep in a 
perturbative region at /i max we can evaluate this factor perturbatively. Two loop calculation is given 
in Ref. Note that scheme dependence is also canceled at this stage and the RGI operator 
becomes scheme independent. 

Our target in this study is the calculation of the first factor Z^J(go, Mmin)- In order to further 
reduce the computational cost we use a relation that Zy = Za implied by the chiral symmetry 
of DWQCD in SF scheme H], together with another one that Zyy^AA — z va+av> which will be 
numerically checked later, and adopt the following definition throughout this paper, 



z v (go) 



2. Renormalization conditions for four-quark operator 

Since we utilize the step scaling function (SSF) obtained by the Alpha collaboration [Q] as 
an intermediate RG running factor from /Vm to /i max , the same renormalization scheme should 
be adopted for our calculation of Z^J(go,Mmin)- The renormalization condition is given by the 
following correlation function [§] 

^r A r fl r c (*o) = ^zi^a^as^Ova+M^^c]), (2-1) 

where subscripts 1 ~ 5 represent quark flavours and Ova+av is the parity odd four quark operator 
which consists of four different flavours from 1 to 4. Boundary operators ffij and ffjj are given in 
terms of boundary fields £ and ^ as 

^■[r] = fl 6 £6(5)r0(y), W =a 6 ££(*)rcj()0. (2.2) 

xy xy 

Due to the SF boundary condition for fermion fields the boundary operator should be parity odd 
and we have two independent choices T = 75 and Y = jk(k= 1,2,3). For the correlation function to 



be totally parity-even we need at least three boundary operators in (2.1). Logarithmic divergences 



in boundary fields £'s can be removed by the boundary -boundary correlation functions; 



/l = -^(^ 2 [75]^2l[75]>, fcl = -^I<^2M^2lM>. (2-3) 

k=l 



We adopt the following three choices [0|, whose perturbative expansions behave reasonably well 

h ±( v x _ ^W,7s]( x o) ,±, \ _ IlLl^s.Tt.Ttl^o) ,±/ \ _ jlLl^[75.7A-:7i]( x o) n 
n \ W ~ 3/2 ' "3 \ x 0) - 3/2 ' n l W - 1/2 ' K - L ^> 

J I J\ Jl M 

We call three renormalization schemes defined through these correlation functions as scheme 1,3, 
7 according to the Alpha collaboration 

We impose the following renormalization condition 

Z$A+AvA8o,H)hf(x = L/2; g0 ) = hf {bee \x = L/2) (2.5) 
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N L 


6 


8 


10 


12 


14 


16 


18 


P 


2.4446 


2.6339 


2.7873 


2.9175 


3.0313 


3.1331 


3.2254 




1.22(2) 


1.32(2) 


1.35(2) 


1.39(2) 


1.40(2) 


1.41(3) 


1.42(3) 



Table 1: j5 values which gives the same box size 2L max for each lattice sizes Nl- 



where s labels the scheme. This means that the renormalized correlation function should coincides 
with that at tree level in the continuum at the middle of the box xq =L/2. The renormalization scale 
at the low energy (hadronic scale) is introduced by the maximum box size 1 / /i m j n = 2L max , where 
L max is defined through renormalized coupling g 2 (l/L max ) = 3.480 in the SF scheme. This box 
size corresponds to L max /ro = 0.749(18) [||] in the continuum limit, so that /Vin = 1 /2L max ~ 263 
MeV using the Sommer scale ro = 0.5 fm. 

3. Numerical simulation details 

The CP-PACS collaboration has calculated the lattice bare Bk in quenched DWQCD with 
Iwasaki gauge action at the domain wall height M = 1.8 and the fifth dimensional length N5 = 16 

In order to renormalize this Bk we adopted the same lattice formulation as in the above. The 
SF formalism for Iwasaki gauge action is given in Ref. [g]. For a domain-wall quark we adopted 
the orbifolding construction [^] to realize the SF boundary condition. We employ the same size for 
the temporal and the spatial directions T = L. We take the mass independent scheme in massless 
limit where all the physical quark masses are set to zero. 

Three lattice spacings j8 = 2.6, 2.9 and 3.2 1 are employed corresponding to aT x ~ 2, 3 and 4 
GeV in the previous simulation [§]. In order to cover these three j8's we take 7 lattice sizes. At 
each lattice size we tune /3 to satisfy oNl = 2L max = 1.498ro using the following fit formula [||] 

^(j-) =-2.193-1.344( j 6-3)+0.191( j 3-3) 2 , (3.1) 

which covers 2.456 < j3 < 3.53. Lattice sizes and corresponding j3 values are listed in Table |T[ 

Quenched gauge configurations are generated by using the HMC algorithm. First 2000 trajec- 
tories were discarded for thermalization. We calculate the correlation functions on each configura- 
tion separated by 200 trajectories. We employed 500 to 1000 configurations in this paper. 

4. Non-perturbative renormalization of B K 

In this section we evaluate the renormalization factor 3fB K ;s(go)-> which convert the lattice bare 
#«"(go) in DWQCD to the RGI Bk- Combining our renormalization factor Zg^. s [go,fi n ^ n ) with 
the RG running factor Z" +AV (oo iJ u max )Zyj , +AV (/x max ,^ m i n ) given by the Alpha collaboration, we 
obtain the renormalization factor ^B K ;s{go) at each /3. A result is given in Table [l] for scheme 1. In 
order to obtain the renormalization factors at j8 = 2.6, 2.9 and 3.2 we fit it in a polynomial form, 
%B K -Ago) =a s + b s (P-3) + c s (fi - 3) 2 . 

'The data at /3 = 3.2 is new and not published in [§] 
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m v a m v a 

Figure 1: Scaling behavior of RGI Bk (left panel) and B/f(MS,2GeV) (right panel). Our results are shown 
by open circle (scheme 1), open up triangle (scheme 3) and open diamond (scheme 7). Previous results with 
tmQCD [^] and DWF [^[ are also shown by open down and right triangles. In the right panel result of the 
KS fermion [|7| and DWF [||| with perturbative renormalizations are also given. The continuum results are 
slightly shifted to avoid overlaps. 



Multiplying it to the bare Bjc(go) we obtain the RGI Bk, whose scaling behavior is shown 
in the left panel of Fig. [jj We also evaluate the renormalized B K in MS scheme with naive di- 
mensional regularization (NDR) at a scale /I = 2 GeV. The scaling behavior of B^ s (NDR, 2GeV) 
is given in the right panel of Fig. ||. Since the scaling violation is negligible, it is reasonable 
to take the continuum limit by a constant extrapolation. We arrival at B K = 0.773(7) ( ) and 
fifS(NDR,2GeV) = 0.557(5) (+j 4 ), where central values are taken from the scheme 1, the first 
parenthesis gives statistical error, and upper(lower) error in the second parenthesis denotes differ- 
ence between schemes 1 and 7(1 and 3). These are main results in this study. 



5. Additional investigations 

We perform two additional calculations. One is to study the scaling behavior of SSF defined 
by Ly A+AV . s (u,a/L) = Z ™+ Av ^ 0,a /^ I j n DWQCD. The other is to investigate a size of 

chiral symmetry breaking effects in renormalization factors. 

5.1 Scaling behavior of SSF at u = 3.480 

We investigate the SSF at strong coupling, u = g 2 (1/L max ) = 3.480, for four different j3's. In 
the left panel of Fig. || we show a dependence of our SSF (open square) and the continuum limit of 
the Alpha collaboration ^ (star). We find that scaling violation of our SSF is large and oscillating 
at M = 1.8. We speculate that this bad scaling behavior is caused by the 6 \a) boundary effect in 
the SF scheme of DWQCD, not by the G{d) bulk chiral symmetry breaking effect. To see this we 
calculate the SSF at tree level, where N5 — > °° limit is already taken. We plot its scaling behavior 
by open circles in the left panel of Fig ||, where the similar sort of oscillating behavior is observed 
at M = 1.8. To exclude a possibility that the large scaling violation is caused by the bulk chiral 
symmetry breaking effect, we calculate ./V5 dependence directly atu = 3.480. Indeed comparisons 
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L (u,a/L) vs. a/L with g (L)=3.480 



□ CP-PACS (M=l. 8) 

▲ CP-PACS (M-1.8 tree level imp.) 

Q CP-PACS (M=1.4) 

yfc continuum (Alpha) 



— m- 
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Z (M=1.4J 
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Figure 2: Scaling behavior of SSF for Zva+av-,i (left panel) and for Z^(right panel). Open squares are 
results for M = 1 .8 and open circles are for M = 1 .4 without improvement. Filled diamonds denote results 
for M = 1.8 with tree level improvement. Star symbol denotes the continuum limit. 



between N$ = 8 and A^5 = 16 for IVa+av^ = 4) and between ./V5 = 32 and ./V5 = 16 for 
LvA+AV;s(u,L/a = 6) show no N5 dependence within statistical errors. 

We find that the scaling behavior of the tree SSF is much improved at M = 0.9, as shown in 
the left panel of Fig [| This suggests that the scaling behavior of the non-perturbative SSF is also 
improved when the renormalized domain-wall height is nearly equal to unity. We first evaluate the 
SSF at M = 1.4, where tadpole improved domain-wall height is nearly equal to unity. Result is 
given by open circles in the left panel of Fig |2| We find good scaling behavior at M = 1.4, and the 
linear continuum extrapolation is consistent with the continuum limit of the Alpha collaboration. 

For another improvement we introduce a tree level improvement, where the continuum tree 
level correlation function hf^ bee \xo) in the renormalization condition (2.5) is replaced to that on 
the lattice hf (^o) LAT - We expect that the large scaling violation of the SSF is partly canceled by 
the tree SSF. We calculate the tree level improved SSF , where the tadpole improved value M ~ 1.5 
of M = 1.8 at each j(3 is used for hf^ ree \x ). The result is given by the filled diamond in the 
left panel of Fig ^[ We find that the scaling behavior is improved, so that we can perform a linear 
extrapolation to the continuum limit, whose value is consistent with that of the Alpha collaboration. 

Furthermore we calculate the SSF of Bv, defined by £R„(w,a/L) = ^fe°' a / 2L ) resu it 1S 

J KK ' ' ' Z BK {g ,a/L) 

plotted in the right panel of Fig ||by open square (M = 1.8) and open circle (M = 1 .4). We find that 
a large scale violation in is partly canceled in ~L Bk at M = 1.8. Linear continuum extrapolations 
using data at finest three lattice spacings give consistent values between M = 1.4 and M = 1.8. 



5.2 Chiral symmetry breaking effect 

We check whether the chiral relation that Zw+aa = Zva+av we assumed is realized or not in 
our DWQCD. If the chiral symmetry were exact we would have a chiral WT identity 

(OvA+AV^[Q) S = {O VV+ AA0[Q)S (5-1) 

under chiral rotation of the first flavour 8q\ = —ijsq, 5£i = ijs^i, SC{ = i7%Ci» Wnere ^[Q is a 
chiral rotation of some boundary operator (pT2|). We get a relation Zw+aa = Zva+av from this 
WT identity. On the other hand, the domain-wall fermion action is not invariant under the chiral 
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TREE LEVEL E. (u,a/L) 
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O ^VA+AV - 




A Z w+AA 
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Figure 3: Left panel is a scaling behavior of tree level SSF. Open circle shows results with M = 1.8 and 
open square shows those with M = 0.9. Right panel is a comparison between ZvA+Av(go, Mmin) (open circle) 
and Zw+AA;i (go, Mmin) (open triangle) as a function of time t . We adopt scheme 1 at the hadronic scale /x m ; n 
on L/a = 16 lattice. 



rotation as Sjwf — > Sdwf + Y, where Y = xj/Xxyis the chiral symmetry breaking term at the middle of 
the fifth dimension. Therefore the WT identity becomes (Ova+av&[C])s = (Ovv+aa0[Q)s+y 
(Ovv+AA^[Q)s, an( l we estimate a possible chiral symmetry breaking effect by directly comparing 
(Ova+av0[Q)s with (Ovv+aa^[Q)s- We evaluate the renormalization factor Z V v+AA(go,Llmm) 
using the chirally rotated correlation function with 100 configurations at each j3. The right panel of 
Fig H shows a time dependence of Zva+av and Zw+aa at L/a = 16 We observe a good agreement 
between them, and a similar results are obtained at other L/a. Therefore the relation Zw+aa = 
Zva+av is realized in our simulations. 

This work is supported in part by Grants-in-Aid for Scientific Research from the Ministry of 
Education, Culture, Sports, Science and Technology(Nos. 18740130). 
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